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Abstract. In this article we study the long-time behaviour of a system of nonlinear 
Partial Differential Equations (PDEs) modelling the motion of incompressible, isothermal 
and conducting modified bipolar fluids in presence of magnetic field. We mainly prove the 
existence of a global attractor denoted by A for the nonlinear semigroup associated to the 
aforementioned systems of nonlinear PDEs. We also show that this nonlinear semigroup 
is uniformly differentiable on A. This fact enables us to go further and prove that the 
attractor A is of finite-dimensional and we give an explicit bounds for its Hausdorff and 
fractal dimensions. 



1. Introduction 

Magnetohydrodynamics (MHD) is a branch of continuum mechanics which studies the 
motion of conducting fluids in the presence of magnetic fields. The system of Partial Dif- 
ferential Equations in MHD are basically obtained through the coupling of the dynamical 
equations of the fluids with the Maxwell's equations which is used to take into account 
the effect of the Lorentz force due to the magnetic field (see for example |9j). They play a 
fundamental role in Astrophysics, Geophysics, Plasma Physics, and in many other areas in 
applied sciences. In many of these, the MHD flow exhibits a turbulent behavior which is 
amongst the very challenging problems in nonlinear science. Due to the folklore fact that 
the Navier-Stokes is an accurate model for the motion of incompressible and turbulence 
in many practical situation, most of scientists have assumed that the fluids in these MHD 
equations follow the Newtonian law in which the reduced stress tensor T(f (u)) is a linear 
function of the strain rate £ (u) = \ ( L + L T ) , L = Vu, where u is the velocity of the 
fluids. In this way we obtain the conventional system for MHD equations which has been 
the object of intensive mathematical research since the pioneering work of Ladyzhenskaya 
and Solonnikov [21]. We only cite [39], [38], [T7], [13] for few relevant examples; the reader 
can consult [18] for a recent and detailed review. 

There are two major problems that arise when studying MHD equations or hydro- 
dynamic equations. First, it is well known that the uniqueness of weak solution of the 
three-dimensional NSE and MHD equations is a very difficult problem. As it is always 
not possible to prove the existence of global attractor in the case of lack uniqueness of 
solution, this is an obstacle for the investigation of the long-time behavior which is very 
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important for the understanding of some physical features (such as the turbulence in hy- 
drodynamics) of the models. We refer, for instance, to [I], [33], and [41 j for some results in 
this direction for the case of autonomous Navier-Stokes and other autonomous equations 
of mathematical physics. Second, there are a lot of conducting materials that cannot be 
characterized by Newtonian fluids. To overcome these two problems one generally has to 
use other model of fluids or some perturbation of the Newtonian law. This has motivated 
scientists to consider (conducting) fluids models that allow T being a nonlinear function 
of S (u). Fluids in the latter class are called Non-Newtonian fluids. We refer for exam- 
ple to the introduction of Biskamp's book [7] for some examples of these Non-Newtonian 
conducting fluids. In [23] and [Mj, Ladyzhenskaya considered a model of nonlinear fluids 
whose reduced stress tensor T(£(u)) satisfies 

T(£(u)) = 2(e + W) |Vun^, r > 0. 

dxj 

Since then this model has been the object of intense mathematical analysis which have 
generated several important results. We refer to [2], [26] for some relevant examples. In 
(14] the authors emphasized important reasons for considering such model. Despite its 
mathematical success the Ladyzhenskaya model has received a lot of negative criticisms 
from physicists. Indeed the Ladyzhenskaya model is a mathematical model used to over- 
come the lack of uniqueness for the NSE; it does not have really a physical meaning as 
it does not satisfy some basic principles of continumm mechanics (the frame indifference 
principle) and thermodynamics (the Clausius-Duhem inequality). Necas, Novotny and 
Silhavy [31], Bellout, Bloom and Necas ]2' has developed the theory of multipolar vis- 
cous fluids which was based on the earlier work of Necas and Silhavy |3Uj . Their theory 
is compatible with the basic principles of thermodynamics such as the Clausius-Duhem 
inequality and the principle of frame indifference, and their results up to date indicate 
that the theory of multipolar fluids may lead to a better understanding of hydrodynamic 
turbulence (see for example [5]). Bipolar fluids whose reduced stress tensor T(f (u))is 
defined by 

T(£(u)) = 2/x (e + |£(u)| 2 )-f £(u) - 2 Ml A£(u), (1) 
form a particular class of multipolar fluids. If < a < 1 then the fluids are said shear 
thinning, and shear thickening when a < 0. 

Let SI be a simply-connected, and bounded domain of R n (n = 2, 3) such that the 
boundary <9ft is of class C°°. This will ensure the existence of a normal vector n at 
each of its point. In this article we are aiming to the long-time behaviour of the system 
of nonlinear PDEs representing the motion of a conducting nonlinear bipolar fluids in 
presence of magnetic fields. This system is basically given by 

' ^ - divT + u- Vu + /ub x curlb + VP = / in ft x (0,oo), 
f£ + (S curl curl b + fj,u ■ Vb - fjh ■ Vu) = in ft x (0,oo) 

div u = div b = in ft x [0, oo), , . 

V \ ( 2 ) 
u = TijirijWi = on <9ft x (0, oo), 

b • n = curlb x n = on <9ft x (0, oo), 

iU (0) = u , b(0)=b inft, 

where u = (uf, i = 1, . . . , n), b = (bf,i = 1, . . . , n) and P are functions defined on ft x [0, T], 
representing, respectively, the fluid velocity, the magnetic field and the modified pressure, 
at each point of ft x [0, oo]. S and \i are positive constants depending on the Reynolds 
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numbers of the fluid and magnetic fields, and the Hartman number. The quantities Uo 
and bo are given initial velocity and magnetic field, respectively. The vector n represent 
the normal to d£l and r%ji is defined by 

d£ij(u) , 

Tiji = q ,i,3,l e {!,••,«}. 

Finally, T designates the extra stress tensor of the Non-Newtonian fluid which is defined by 
(HJ). We suppose throughout that s, [Ao,(J,i are positive constants and a G [0, 1). Hereafter 
we set 

r(u) = / u {e+\£(u)\ 2 y f ■ 

The structure of the nonlinearity of problem ([2]) makes it as interesting as any nonlinear 
evolution equations of mathematical physics such as the conventional MHD or the Navier- 
Stokes equations. 

When b = then ([2]) reduces to the PDEs describing the motion of isothermal incom- 
pressible nonlinear bipolar fluids which has been thoroughly investigated during the last 
two decades. Existence and uniqueness results of weak solution for nonlinear bipolar fluids 
with homogeneous boundary condition were given in [3], [3]. Existence of unique solution 
and asymptotic stability of the solutions for the case for nonhomogeneous boundary con- 
dition were investigated in [6]. The long-time behaviour of (J2J) with b = is investigated 
in [5] for a £ [0, 1) and in [8] for a € (— oo, 0). Both of the authors of [5] and [8] proved 
the existence of a global attractor; they also show that the flow is finite-dimensional by 
giving explicit bounds for the Hausdorff and fractal dimensions of the global attractor. 
Existence results for multipolar fluids are also given in [27] and |31j . These are just ex- 
amples of relevant work related to ([2]) with b = 0, the reader can consult [29] and [16] for 
detailed and recent reviews of results for nonlinear bipolar fluids. 

For a = 0, fi\ = 0, ([2]) is reduced to MHD equations which has been the object of in- 
tensive mathematical research since the pioneering work of Ladyzhenskaya and Solonnikov 
|21j . We only cite [3IJ]> |38j . |17j . |13| for few relevant examples; the reader can consult 
|18] for a recent and detailed review. Assuming that fj,± = Samokhin studied the MHD 
equations arising from the coupling of the Ladyzhenskaya model with the Maxwell equa- 
tions in [37] , [3S] > [M] j and [3S] . In these papers he proved the existence of weak solution 
of the model for a < 1 — ^rj- Later on Gunzburger and his collaborators (see [19] and 
|20j) generalized the settings of Samokhin by taking a fluid with a stress tensor having a 
more general structure. The authors of the last two papers analyzed the well-posedness 
and the control of ([2} still in the case where fio = and a < 1 — 

To the best of our knowledge the preprint [32] is the only work treating the case {/ii ^ 
0,b^0,ae[l-^, 1)}. It was proved in [32] that if (u ; b ) G I and / £ V*, then © 
with e = 1 (but the author's result is still true for any e > 0) has at least a global weak 
solution. This means that there exists a couple (u; b) such that 

. (u;b) eL~(0,oo;V)nL2 oc (0,oo;M), 
• (u; b) satisfies 

*) - - ff) + < r <»e» £ M. ^» + (-W • **>.*> 

-^(b(t)-Vb(t),0) = (/,0), (3) 
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and 



db(t) 




S (Ab(i), V) + M(u(t) • Vb(t), VO - A*(b(t) • Vu(t), V) = 0, 



(4) 



for any couple = (<f>; ip) of smooth solenoidal functions and t G [0, T], T > 0. In 
([3|) and throughout this work summations over repeated indices are enforced. 

The author also analyzed the long-time behaviour of solutions of ([2]) with a € [1— ^Bj; 1)- 
He mainly proved the existence of trajectory attractor for the translation semigroup acting 
on the trajectories of the set of weak solutions and that of external forces. In the case 
when a E [1 — ^pg, 1) we are not sure of the uniqueness of the weak solution, therefore 
the author of [32] follows the method in [10] (see also, [11] and [12]). 

For < a < 1 (which correspond to 1 < p < 2 in the paper [32] ) it is not difficult 
to check that the weak solution is unique. We will formally check this statement later 
on. This fact allows to define the nonlinear semigroup {§(i);i > 0} of the solution of ([2]). 
It is our purpose to study the long time behavior of this semigroup. The semigroup is 
formally defined by " let t > to each initial value (uo; bo) € EI we associate an element 
§(i)(uo; bo) = (u(t); b(i)) of EI which is the unique weak wolution of ([2])." Our results are 
as follow: 

(1) We give the existence of a compact global attractor A for S(t). For this purpose 
we need to find absorbing sets in HI and in V. Since V is compact in EI then we 
can define from classical argument the existence of the compact global attractor 



(2) We show that the semigroup S(i) is differentiable with respect to the initial data 
(uo; bo) € A. This result is crucial for establishing the next result. 

(3) We give an estimate for the bounds of the Hausdorff and fractal dimensions of A. 

To establish these results we mainly follow the idea in [5] and the classical results for the 
investigation of long-time behaviour of dissipative PDEs presented in jl] , [33J and [H] for 
example. We should mention from the very beginning that every calculation we performed 
is formal, but we can check them rigourously by using the Galerkin approximation and 
pass to the limit as it was done in [32]. Also, even if we drew our inspiration from [1], [5], 
[15] . [33j, [38] and [H] the problem we treated here does not fall in the framework of these 
main references. Besides the usual nonlinear terms of the conventional MHD equations it 
contains another nonlinear term which exhibits the non-linear relationships between the 
reduced stress and the rate of strain £ (u) of the conducting fluids. Because of this, the 
analysis of the behavior of the MHD model ([2]) tends to be much more complicated and 
subtle than that of the Newtonian MHD equations and the model in [5]. Hence, we have 
had to invest much effort to prove many important results which do not follow from the 
analysis in the aforementioned papers. 

The organization of this article is as follows. We introduce the necessary notations for 
the mathematical theory of ([2]) in the next section. In Section 3 we prove the existence of 
absorbing sets in EI and in Y which enables us to show the existence of a global attractor 
for the nonlinear semigroup of solutions to ([2]). The uniformly differentiability of §(i) on 
A is establisehd in Section 4. In the last section we give explicit bounds for the fractal 
dimension of the global attractor A. 



A. 
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2. Preliminary 

We introduce some notations and background following the mathematical theory of 
hydrodynamics (see for instance [10]). For any p G [1, oo), L p (£7) and W m ' p (0) are the 
spaces of functions taking values in R n such that each component belongs to the Lebesgue 
space L p {n) and the Sobolev spaces W m,p (n), respectively. For p = 2 we use HI m (f2) 
to describe W rn,p (Q). The symbols | • | and (., .) are the L 2 -norm and L 2 -inner product, 
respectively. The norm of W 2,m (f2) (resp., W p,m ) is denoted by || • || m (resp., ||-||p, m ). As 
usual, C§°(0) is the space of infinitely differentiable functions having compact support 
contained in n. The space Wg' m is the closure of C§°(0) in W p,m . Now we introduce the 
following spaces 

Vi = {ueCg°(O):divu = 0}, 

Mi = {uGL 2 (ll) :divu = 0,u-n = Oondft}, 

V 1 = Jug U 2 (Q) : divu = 0,u = ^ = on dtt\ . 
[ an J 



We also set 



V 2 = {b G Cg°(ft) : divb = 0;b • n = on dn} 



H2 



the closure of V2 in L 2 (S1), 



V 2 = {b € H 1 ^) : divb = 0;b • n = on dn} , 

V 3 = {b G M 2 (n) : divb = 0;b • n = curlb x n = on dn} 



Note that 



111 = M12- 

112/ 



We equip the space Vi with the norm 1 1 - 1 1 2 generated by the usual H (O)-scalar product. 
On V2 we define the scalar product 

((u, v))v 2 = (curl u, curl v), 

which coincides with the usual scalar product of H 1 ($7). Hence, from now on ||b||i we will 
denote the norm of b £ V2 in H 1 (r2) and ||b||v 2 = |curlb| as well. 
Let 

V = Vi x V 2 , 
H = Hi x H 2 . 

The space EI has the structure of a Hilbert space when equipped with the scalar product 

(*,¥) = (u,v) + (b,C), (5) 

for $ = (u;b),\P = (v;C) G H. 

The space V is a Banach space with norm 

||$||v = ||u|| 2 + ||b||i, (6) 

for $ = (u; b) G V. 

For any Banach space X we denote by X* its dual space and (</>, u) the value of <f> G X* 
on u G X. 

For any any u G W k+1,p , k > 0, we set 

£(u) = I T(Vu) + (Vu) Tl 



(In 


d k £{u) 




dxp 1 .dxp l 



P dx) , (7) 
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Let us recall the following results whose proofs can be found in |28j . 

Lemma 2.1 (Korn's inequalities). Let 1 < p < oo, (3 = (/3i, . . . , /?;) be a multi-index 
such that k = Ya=iPi> an d ^ C K n be a bouned domain with smooth boundary. Then 
there exists a positive constant K(£l) such that 

K(n)\\u\\ p , k+1 < 

for any u G W fc+1,p . 

The following Poincare's inequality is very crucial 

Ai|u| 2 < [|u||f, (8) 

for any u G H 1 (J7). This inequality holds with 1 1 u. 1 1 2 for u G EI 2 (f2). 

We also recall that there exist two positive constants Ki(£l), /^(O) such that 

Xi(0)||b||^ 2 < |Ab| 2 < K 2 (n)\\b\\ m 2, (9) 

for any b G V 2 nrfn). 
There also holds 

|V0| 2 < H 2 + |curl(/>| 2 , (10) 
for any divergence free function eft satisfying ((f) ■ n)\ gn = 0. We refer to [22] or [25] for the 
proofs of © and dlHl). 

3. Existence of the Global attractor A 

In this section we will investigate the existence of a set A C HI which is the global 
attractor of 8(i) in H. But first of all we recall some definitions which are taken from |41j . 

Definition 3.1. A set A C EI is called an attractor for the semigroup §(i) if it enjoys the 
following properties: 

• A is invariant; that is, §(i).4 = A, for any t > 0, 

• A possesses am open neighborhood U such that for any uo G U 

lim inf |§(i)uo — y\ = 0. 

t-S-oo y£A 

The set A C EI is said to be a global attractor if it is compact attractor and attracts any 
bounded sets of H. 

The notion of an absorbing set is very important for the investigation of the existenceof 
a global attractor A. We also recall the definition of an absorbing set B. 

Definition 3.2. Let B be a bounded set of EI and U be an open set containing B. We say 
that B is an abosrbing in U if for any bounded set Bq C U, there exists a time to = ^o(£>o) 
such that S(t)Bo C B, for any t > to- 
After these definitions we formulate our first main result. 

Theorem 3.3. The semigroup S(t) : EI — > EI, t > has a global attractor A which is 
maximal wrt to the inclusion in EI. 

To prove this theorem we start by finding absorbing sets in EI and V. 

Lemma 3.4. There exists a ball C EI which absorbs any bounded set ofW. 
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Proof. The set V is dense in H, so it is permissible to take <f> = u in Q. From this 
operation we obtain that 

\j t ^ 2 + "> (^^r 1 ' ^r) + (r(uw)£(u (t )),£(u(0)) 

- M (b(*) • Vb(t), u(t)) - (/, u(t)) = 0, (11) 

where we have used the well-known fact that (w • Vv, v) = for any w € H l (fl) and any 
divergence free function v £ i^ 1 (S7) . Using a similar reasoning we also have that 

^|b(t)| 2 - S(Ab(t),b(t)) - M (b(t) • Vu(t),b(t)) = 0. (12) 
Now summing up ([TT]) and (fT2|) side by side yields 

-5(Ab(i),b(i))-(/,u(t)) = 0. (13) 
Here the property (b • Vb, u) = — (b • Vu, b) was used. Since 

(r(u)£(u),£(u)) >0, 

and 

-5(Ab,b) > 5||b||f, 

it follows from (fl~3l) that 



~ (|u(t)| 2 + |b(i)| 2 ) + Mr(n)||u(*)||! + S\\h{t)\\\ < \g\\u{t)\. (14) 

Letting vq = mm(fj,iK(Q), S) we can see from (fl~4"l) along with Poincare's inequality ([HJ 
and Young's inequality that 

(|u(i)| 2 + |b(t)| 2 ) + ^ Ai (|u(t)| 2 + |b(t)| 2 ) < 2<T 1 | 5 | 2 + 52- 1 (|u(t)| 2 + |b(t)| 2 ) , 

for any 5 > 0. Choosing <5 = vq we derive from the last estimate that 

(|u(t)| 2 + |b(t)| 2 ) + ^oAi (|u(t)| 2 + |b(i)| 2 ) < -\g\\ (15) 

OA, U(j 

Setting ^i = 7^ and y(t) = \u(t)\ 2 + |b(t)| 2 with y(0) = |u | 2 + |b | 2 , the estimate ([TSJ is 
equivalent to 

+ v y(t) <^i|/| 2 , 



(it 

which together with Gronwall's lemma implies that 

y(t)<y(0)e-^ + ^S-(l-e- nt ) 

<l/(0)e-^* + ^L. 
Since e _Vl * — >■ as f — >■ oo, we can find a time to = *o( u o> bo) such that 



y(t) < 2 -^l, (16) 
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for any t > to- This means that the ball in M 

= {zei; fz| 2 < p\}, 

- J is an absorbing set in H. □ 

Having found B^ 1 we prove that there is also a ball B^ C V which attracts all bounded 
sets of M. We postpone the proof for the next lemma. For now we establish an additional 
estimate that is going to play an important role in the aforementioned claim. From (|14|) 
we see that 

^ + 2u (||u(t)||I + ||b(t)||?) < 21/1111(01, (17) 
for any t > 0. Integrating (|17p over [t,t + r], with r > arbitrary, yields 

y(t + r) + 2i, (||u( 8 )||l + ||b(a)||?) ds < y(t) + 2|/| ^ |u( S )|ds. 
Whenever t >to, we have that 

y(t + r) + 2^ jT^ (||u( S )||| + ||b(s)|| 2 ) ds < p\ + 2|/| Pl r. (18) 
Out of (|18p . we easily derive that 

t+r (||u( S )||2 + ||b( S )||?)d S < Ko (r), (19) 



where 

«o(r) = 2^(p? + 2|/|pir). 

We will now show that there is also a ball By C V which absorbs all bounded set of HI. 

Lemma 3.5. The semigroup §(t) has a bounded absorbing set By C V. 

Proof. The proof of the result will consist of several steps. 
Step 1: Uniform estimate wrt t of b in V2 
Taking ip = — Ab(t) in Q we obtain that 

- (^-,Ah(t)\ + S\Ah(t)\ 2 = M [(u(t) • Vb(t), Ab(i)) - (b(i) • Vu(i), Ab(t))] , 
which is equivalent to 

™| curlb(t)| 2 + S\ Ab(t)| 2 = /x [(u(t) • Vb(t), Ab(t)) - (b(t) • Vu(t), Ab(t))] . (20) 

Since | curlb(i)| 2 is the same as ||b(i)|| 2 for b(i) € Vi, we can rewrite (|20p in the following 
form 

^l|b(t)|| 2 + S|Ab(t)| 2 = fx [(u(t) • Vb(i), Ab(t)) - (b(t) • Vu(t), Ab(i))] . (21) 



Now we want to estimate the right hand side of (|2ip . To do so we mainly use that fact 
that if u e Vi, then u G L q (Q) for any 2 < q < 00 and Vu E L s (ti) with 2 < s < ^ 
(2 < s < 00 if n = 2.) For the first term we have 

fi(u • Vb, b) < ^ci|u|loo |Vb||Ab| 

< 5" 1 (/xci) 2 ||u||l|Vb| 2 + <5|Ab| 2 . ^ 



(24) 
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Owing to CLU and ([22]) we have that 

/i(u-Vb,Ab) < (T^^^llullillbH? + 5\Ab\ 2 . (23) 
By Holder's inequality we see that the second term verifies 

/x(b • Vu, Ab) < /ic 3 |b| L 4|Vu| L 4|Ab| 

< ^ 1 (^c 4 ) 2 ||b||?||u|^ + 5\Ab\ 2 . 
Let us set v\ = 2<5~ 1 / u min(c4, C2). Owing to ([2"3"]) and ([2~4"]) we infer from ([2"T]) that 

^l|b(t)||? + 5|Ab(t)| 2 < v 2 \\h{t)\\l\\u{t)\\l + 25\Ah{t)\ 2 , 

which with an appropriate choice of <5, let us say 5 = j , implies that 

lllbWiif + ^AbWi^^iiMt)!! 2 !^)!! 2 . (25) 

Because of (|19p we can use the uniform Gronwall lemma (see, for instance, [41]) in (j25|) 
to get an uniform estimate wrt t of b in V2 which does not explode too fast in time. This 
estimate reads 

|| b(t + r) ||2<fWlO e ^O (r) 

r 

for any r > and t > to- This implies that for any r > there exists Ki(r) = 
ifi^W e <Wr) > o such that 

l|b(i)|| 2 <Ki(r), (26) 

for any t > to + r. 

Step 2: Uniform estimate wrt t of Ab in L 2 (t, t + r; L 2 ({1)) 
For t > to + r we obtain from (|25p and (|26p that 

»t+r r t+r 



\\b(t + r)\\ 2 + S [ r \Ah(s)\ 2 ds < v lKl (r) [ " ||u(s)|| 2 ds + ||b 
Jt Jt 

which along with (fl~9|) and ([26]) enables us to deduce that 

ft+r 

S J \Ab(s)\ 2 ds < Kl (uiK {r) + 1) . 

We deduce easily from this last estimate that for any r > 

+ \Ah{s)\ 2 ds < K 2 (r), (27) 

for any t > to + r. This inequality will be very useful for the next step. 
Step 3: Uniform estimate wrt t of u in Vi 

We want to establish an uniform estimate of du/dt in HI. To shorten notation we will set 
u t = du/dt. We take 4> = u t hi © and we get 

|2 , fd£(u(t)) dS(n t (t))\ 



|ui(t)I + 111 \ dx k ' dx k ) + ( r ( u )^( u (*))^( u *)) + ( U W • Vu W> U *W) (28) 

= Kb(t).Vb(f),u ( (i)) + (/,u t (t)). 

Let E be the potential defined by 

£(|£| 2 )= / ^o{e + s)- a/2 ds. 
Jo 
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it is not difficult to see that 

|s(|^(u(t))| 2 ) = Mo (e + \E (u(i)))- Q / 2 £ (u(t))£ (u t (i)), 

and 

'dS(u(t)) d£(u t (t))\ d fd£(u(t)) d£(u(t)) 



\u t (t)\ 2 + j t ^J^(£(u(t))dx + ^ 



dxk dxk J dt \ dxk dxk 

Hence we deduce from (1281) that 

d£(u(t)) d£(u(t)) 

dx k ' dx k J J (29) 
< -(u(i) • Vu(t), u t (t)) + (i(b(t) ■ Vb(t), u t (t)) + \f\\u t (t)\, 

which immediately implies that 

< -(u(i) • Vu(i), u t (t)) + M (b(t) • Vb(t), u t (t)) + i|/| 2 , (30) 

</ 1 + / 2 + I|/| 2 . 



(31) 



For ii we have 

|/i| < C5|u(t)U<x,|Vu(*)||ut(t)|, 

<c 6 r 1 ||u(t)||||Vu(t)| 2 + «5|u t (i)| ; 

for any 5 > 0. We can also check that 

\h\ < ^crlbWi^Vbit)] 2 + 5\u t (t)\ 2 , 

for any 5 > 0. Thanks to Q ||b||| is equivalent to |Ab| 2 for any b £ V2 fl V3. Therefore, 
we can derive from the last estimate and (jlUj) that 

|/ 2 | < 5- 1 c 8 \Ab(t)\ 2 \\b(t)\\j+5\u t (t)\ 2 . (32) 

By choosing 5 = 1/8, we can deduce from (|3T)j) - (|3"2j) that 

l| Ut ( t )|2 + ^) < 8c 6 ||u(i)||i|Vu| 2 + 8c 8 |Ab(t)| 2 ||b(t)|| 2 , (33) 
where we have set 

z(t) = E(£(u(t))dx + m ( dS ^\ d£(u(t))dx k ^j . 

Noticing that 

'dS(u(t)) d£(u(t)) 



and dropping out the term l/4|u 4 (i)| 2 we deduce from ([33]) that 

^<^) + 8c 8 jAb(t)| 2 ||b(t)|| 2 . (34) 
at Hi 
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From here we want to use the Uniform Gronwall lemma, so we need to check that for 
certain t' and any r > there exist positive constants ai, 02, 03 such that 



8£6 

Ml Jt 



/t+r 
\Vu(s)\ 2 ds < ai, (35) 

. |Ab(s)| 2 rfs||b(s)||?rf s < « 2 , (36) 

rt+r 

z(s)ds < as, (37) 



for any t > t' . Thanks to previous estimate we take t' = to + r and infer from (|19j) that 

8c 6 f t+r 2 8c 9 2 

— / Vufs) as < — / u(s) Unas 
Mi Jt Mi Jt 

^ 8C9 

Mi 



, t+r)+2mr( «,M )&+2 ^ ar(u)l , (u(s))|2)rfs 

+25^ " \\u(s)\\ 2 1 ds<2p l \f\+y(t), 



Setting ai = ^-Ko(r), we get (|3"5]) . Invoking ([2U]) and ([2"T|) we see that 

rt+r rt+r 

8c 8 / |Ab(s)| 2 ||b(s)||frfs < 8c 8 sup ||b(a)||f / \Ab(s)\ 2 ds, 

Jt s>t +r Jt 

< 8c 8 Ki(r)K 2 (r). 

Letting ai = 8csKi(r)n2(r), we obtain (|36|) . Now let us deal with ([37]) . Thanks to ([13 
and (1161) we have that 



(38) 



for any t > to + r - Keeping only the second term of the left hand side of (|38p and using 
([T6|) we see that 

To estimate the term involving £(£(u)). For this purpose we notice that the function 
g(s) = po(e + s)~ a / 2 is decreasing for s G [0, 00). Hence 

£(£(u))< /" (U)I (supg(s))ds = ^\£(u)\ 2 . 
Jo s>o 

This implies that for any t > to + 2 

^ (7 E(£(u(*)))) * < ^ |^(u( S ))| 2 d S . 
Invoking korn's inequality and (|19j) we have that 

E(£(u(*)))) d S < J^K(n) K0 (r)r. 
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So putting 03 = pi(\f\ + pi) + j£j2K(Q)no(r)r, we get ([37|) . Now we can apply the Uniform 
Gronwall lemma to (|34p and we get 

z(t + r) < (— + a 2 ) exp(ai), (40) 



r 

for any t > to + r - Korn's inequality along with ([40 p implies that for any 

^ 2 ^jdw)^ +a ^ eMaih <41> 

for any t > to + 2r. 
Let us set 

1 /o 3 



K3(r) = Kl(r) + ^(oj 17" + a v exp(ai) - 

So combining (|26p and (|4ip we obtain that 

l|u(t)||2 + ||b(t)||?<«3(r), (42) 
for any t > to + 2r. The equation (|42p implies that the set 

^ 2 ={z€V: ||z||^<pl}, 



with p2 = \J t3(r), is an attracting set in V = Vi x V2. Here we have set 

II Il2 11 ||2 1 iiull2 
ll z llv = Il u ll2 + ll b lll> 

for any z = (u; b) £ V. □ 

Proof of Theorem \3.3[ It is now standard to prove the existence of the global attractor. 
The set B^ is absorbing in HI and there exists to such that E(t)B^ C -By 2 , for any t > to- 
Moreover, By is compact in HI, hence L)t>t E(t)B^ is compact and we can deduce from 

[4"T1 Theorem 1. 1.1] that A = uj(B^) (= r\t Q >oUt>t ^(t) ) is a compact attractor 
which is maximal wrt the inclusion in HI. This completes the proof of the theorem. □ 

Our next concern is to find an estimate of the bounds for du(A) and df(A), the Haus- 
dorff and fractal dimension of A. This procedure needs that S(t) is Frechet differentiable 
wrt to the initial data. Will show this fact in the next section. 

4. Uniform differentiability of §(t) 

Our second result is about the regularity of the semigroup {S(t) : t > 0} with respect 
to the initial data of the problem ([2]). We recall the following definition which is borrowed 
from 



Definition 4.1. We say that S(t) is uniformly differentiable on A if for every (u; b) G A 
there exists a linear operator £(u; b), such that, for all t > 0; 

|S(t)(v; c) - S(t)(u; b) - £(u; b)(v - u; c - b)| 

SUp : : : —. > 0, 

{(v;c), (u;b)e^: |v-u| + |c-b|<x} |V — U| + |C — b 

as x — > 0, and 

[£(u;b)(v;c)| 

sup sup : — : j— : < 00, for each t > 0. 

(u;b)G-4(v;c):|v| + |c|^0 l V l + l C l 

The main result of this section is stated in the following 
Theorem 4.2. The semigroup S(t),t > is uniformly differentiable on A. 
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Proof. We start the proof by showing that §(i) is Lipschitz continuous. For doing so we 
consider two weak solutions (u; b) and (v; be) associated to the same forcing / and the 
initial data (uo;bo) and (vo;co), respectively. The functions w = v — u and m = c — b 
satisfy 



d£(w(t)) d£(<j>)\ 



dx k dx k J 
+ (v(t) • Vw(t) + w(i) • u(t),<f>) + pi 



+ (r(v(t))f(v(t))-r(u(t))f(u(t)),£(^) 

(b(t)-V^m(t)) + (m(t)-V0,c(t)) 



0, (43) 



and 



(m t ,^)-S(A,^)+fi[ (u(t)-Vm(t)-m(t)-Vu(t)-c(t)-Vw(t)+w(t)-Vc(t),^) ) = 0, 

(44) 



where we have set wj = — — and ni( = — — . Taking 

at at 

equations yields 



w and tp = m in the above 



i|w«)P+»f^i.^a| + l'-.v«<.»-.v«,,.-rm(- ) ) f -«»,.).. f .w„.) 



(w(t).Vw(i),u(t))+/i 



(b(t) • Vw(t),m(f) + m(t) • Vw(t),c(f)) 



~\m(t)\ 2 + 5||b||f = n[ (m(t) • Vu(t),m(*)) + (b(t) • Vw(i),m(t)) 



and 
1 

2 eft 



Let $(i) = |w(i)| + jm(i)| 2 . Summing up side by side and applying Korn's inequality we 
obtain that 



+ (w(t) • Vm(t),c(t)) . 



ld_ 

2di 



$(t) + (||w(i)||l + ||m(t)||i) < |(w(t) • Vw(i),u(t))| + /x 



(m(i)-Vw(i),c(t)) 



(m(t) • Vu(i),m(t)) + (w(t) • Vm(t),c(t)) 



(45) 



To estimate the right hand side of (|45p we will mainly use the fact that if (uo;bo) and 
(vo;co) are an element of A then (u(£);b(t)) and (v(t);c(t)) are element of -By 2 for any 
t > 0. For the first term we have that 

|(w(t).Vw(t),u(t))| <ci||u(t)|||w(t)||Vw(t)|, 

where we have used Holder's inequality and the embedding H 2 (£l) C L°°(Cl). By a similar 
argument we have 

|(m(t) • V,c(t)) + (w(i) • Vm(t),c(t))| < c 3 |m(t)| |Vw(t)| L 4 |c(t)| L 4 

+ ci|m(t)||Vc(t)||w(t)| L oo, 
<Ps|m(t)|||w(0[| 2 ||c(*)||i. 
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We can also check that 

|(m(t)- Vu(t),m(t))| <C3|m(*)||Vu(t)| L *|m(*)| L 4 

< c 3 \ m (t)\\\Mt)h\Ht)h- 

Since (u(i),b(i)) and (v(t),c(t) are element of By , we can find a positive constant cio 
such that the right hand side of ([45 p is bounded from above by 

2/52C10 (|w(t)|||w(t)|| 2 + |m(i)|||m(i)||i + |m(t)||Am(t)|) , 

where we used the fact that |jm||2 and |Am| for m G V2 fl H 2 (£l). Hence 

+ v [\\w(t)\\ 2 2 + ||m(t)|||) < 4p2Ciof (|m(t)| + |w(t)|) (||w(t)|| 2 + ||m 
from which we can infer that 

|*(t) + 2v (||w(t)i|I + ||m(t)||i) < 2x- 1 (8 P2 c 10 ) 2 1>(t) + x(||w(t)|| 2 + ||m(t)|| 2 ), 
for any x > 0. By choosing x = Uq, we deduce from the last inequality that 
-$(£) + 1/0 (||w(t)|[! + ||m(t)||l) < 2x- 1 (8p 2 c w )H(t). 

We put 771 = 2x _1 (8/?2Cio) 2 and use Gronwall's lemma to deduce that 

< (|w | 2 + |m | 2 ) expfat), (46) 

for any i > 0. This show that is Lipschitz continuous wrt to the initial data, which 
show also the uniqueness of the solution of fl2J). 

Let (u;b) a weak solutions of Pj). Using the same argument as in [5] it can be shown 
that the linearization of (J2J) about (u; b) is given by the following system of linear PDEs: 

(^T>A + (r(u(t))f (u(i)) - aA p (u(t))5(u(t)), f (0)) + Ml (£ (u(t)), 5 (0)) 



= m( (M(t) • Vb(i), <p) + (b(i) • VM(t)»l + (U(t) • V0, u(t)) 

+ (u(t)-V0,U(t)), (47) 

" ^ AM W'^ = ~^(( U (*) • Vb(t),V) + (u(t) • VM(t)i) 

- (M(t) • Vu(t), - (b(t) • VU(t), V) J • (48) 

Here the tensor A^-jy is defined by 

A ijkl (u) = (jlq (e + |f(u)| 2 ) 1+a/2 %(u)^(u), (49) 

and it satisfies 

r(u)|5 (U)| 2 - aAi ifc ,(u)%(U)£fc,(U) )daj 

(50) 



>9, f |g(U)l 2 ,,,, \ / ^(U)%(U)^(u)g fc; (u) 

> 2eafi / — — jrrdx + 2(1 - a)uo / — jt, 

Jn (e + \£(M)\ 2 ) l+a/2 ' Jn (e + \£ (u)|T /2 



QUALITATIVE PROPERTIES OF THE MHD EQUATIOSN FOR BIPOLAR FLUIDS 15 

for all £,no > 0, < a < 1. Note that by standard Galerkin method we can show that 
(|47p ~(|48 p has a unique solution (U;M) € L^ c (0,oo;H) n Lf oc (0, oo; V) which implies the 
second condition in Definition 14.11 Now let = v — u — U and ^ = b — c — M. After 
some algebra 



+ a (A ijkl (u(t)) £ (U(t)), £ {#)) + (i f(b(i) ■ V<j>, vP(i)) + (*(t) ■ V<p, b(t)) 
+ (v(t) • V0(t) + 6(t) • Vv(t), 0) - (w(t) • Vw(t), 0) + (m(t) • Vm(t), 0) = 0, (51) 

and 

(^U) - 5 ( A *(*)^) + /*(e(t) • vc(i) - b(i) • ve(t) + v(t) • v*(i) 



- *(t) • Vv(t) + m(t) ■ Vw(t) - U(t) • Vm(t), V> 1 = 0. (52) 
Taking <p = (resp., tp = $>) in the ([51]) (resp., ([52]) ) yields 



r(u(t))f(u(t)),f(e(t))J +a(A< iH (u(t))£(u(t)),£(e(t))) + (e.vv(t),e(t)) 
= ( w (i) • vw(t), e(t)) - //(b(t) • ve(t),f(t)) + • ve(t), b(t)) 

+ (m(t) • Vm(i),9(t)) 
and 



+ (m(i) • Vw(t) - U(t) • Vm(i), *(*))! = (53) 
By using some results from Page 150-151] we see that 

l^dT- + ^( fi )ll (*)lli + ( W • vv(i), e(i)) - (w(t) • vw(t), ©(t)) 

< -/i(W) • V6(i), + (tt(t) • V9(t), b(i)) + (m(t) • m(t), 9(t)) 

(e(t))£ kl (w(t))£ mn {w{t))dx, (54) 

where 

-1 /•! 



ijkimn= / o r ~ f p — (£(u(i) + CTTw(i))) rdrdcr (55) 

J0 JO OOijOOklO mn 
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dt 



+ S\\*{t)\\i<» 



(Q(t) • Vc(t) - b(t) • V0(t), *(*)) - (tf (t) • Vv(i), *(*)) 



(56) 



+ (m(t) • Vw(t) - U(i) • Vm(t), *(t)) 
Now we want to estimate each term in the right hand side of (|54p . We have that 

e(t) ■ vu(t),e(oj < d|e(t)||ve(t)||u(t)| LO o, 
<c(fi)|G(i)|||e(i)|| 2 ||u(f)|| 2 . 

Owing to Young's inequality and the fact that u(t) £ By for any t > and , we obtain 
that 



e(t) • vu(t), e(i) < c(n)V2^~ !©(*)! + *H6(*)ll2 



(57) 



for any x > 0. 

Holder's inequality implies that 



b(t) • ve(t), *(t) + *(t) • ve(t), b(i) 



<2C(n)|b(t)| L 4|ve(t)| L 4|*(t)|. (58) 



Using fact that iT C L 4 (f2) for n = 2,3, and the HV'll.ff 1 an d | curl -0 1 for any tp € V 2 
we deduce from the last inequality that 



b(t)- ve(t),*(t) + *(t).ve(t),b(*) 



<2C(J2)||b(t)|| 1 ||e(t)|| a [*(t)|. (59) 



Since (u; b) £ J3y , we easily see that ||b(i)||i < p 2 - Therefore 



b(*).ve(t),*(t) + *(*)• ve(i),b(t) 



It is not difficult to show that 



< x- 1 (2/xC(n)p 2 ) 2 |#(i)| 2 + x\\e(t)\\ 2 2 . 

(60) 



b(t)-Vm(i),9(t) 



</iC(n)||m(t)||i||e(t)|| 2 |m(*)|. 



Since (w;m) E By , we have |[m||i < p%. Thus, 



b(t) • Vm(t),0(t) 
We also check that 

w(t)-vw(t),e(t) 



< x- 1 (fi P2 C(n)) 2 \m(t)\ 2 + x||G(t)|||. 



(61) 



<c(n)[w(t)|||e(t)|| 2 ||w(t)|| 2 , 

<^ 1 (C(n) P2 ) 2 \w(t)\l + K\\Q(t)\\l (62) 
As far as the term involving ^ijkimn is concerned we have that 

EijMmnSij (©(*)) ffci (w(t)) £ m „ (w(t)) < x- 1 (C(fi)p 2 ) 2 || w(t) ||| + x||9(t) g, (63) 
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where we have used again the fact that 1 1 w(t) 1 1 2 < pi- Putting ([5"7]) - (f63j) in (|54"|) yields 
ld\e(t)\ 2 



+ ^K(n)\\@(t)\\ 2 < K- 1 (c(n) P2 ) 2 \e(t)\ 2 + ^ 2 \^(t)\ 2 + p 2 \ m (t)\ 



+ |w(i)| 2 + ||w(t)|| 2 +5x||GW|| 2 (64) 



dt 



In the next few lines we estimate each term in the right hand side of (|56p . Owing to the 
same argument we used for (|59p and the fact that |Vb| < |curlb| + |b| we obtain that 



G(i) • Vc(t) -b(t) • VG(i),#(t) 



<c(n) |vc(t)|||e(t)|| 2 |*(t)| 



|b(t)||i||e(t)|| 2 |*(t) 



<c(0) (||c(t)||i + |c(0|)||e(t)|| 2 |*(t)| 



+ ||b(t)||i||e(t)|| 2 |*(t)| . 



Since (u; b) and (v; c) belong to -By 2 we have that 



9(t) • Vc(t) - b(i) • V0(t), *(t) 



<3//p 2 C(fi)|*(i)|||G(i)|| 2 

< x' 1 (3fip2C(n)) 2 \V(t)\ 2 + x||Q(i)|||. (65) 
Using the similar idea as used for (|59p and (|65[) we see that 



P 



(*(*)■ Vu(t),*(t)) < pC(n)\*(t)\ Li \Vu(t)\ L *m)\, 

< M c(n)||*(t)||i||u(t)|| 2 |*(t)|, 
< 7 - 1 (mC(0)) 2 |vi/(t)| 2 + 7 ||M/(t)||2. 

Before we proceed further we recall that U = w — G. Hence 

m(t) ■ Vw(i) -U(i) • Vm(t) , = f m(i) • Vw(i) - w(t) • Vm(i) + 9(t) • Vm(t), *(t) 

As for (1661) we can check that 



(66) 



m(i) • Vw(i),f(t) 



< 7 -Vp2<?(n))V(*)| 2 + 7ll*(*)lli- 



(67) 



We also have that 

w(t) • Vm(t),$(t) 



< C(fi)|w(t)|r« |V*(t)||m(t)|, 



< C7(«) ^|[w(0 || a ||i I + ||w(t)|| 2 |m(t)||*(t)|J . 

If (uo;bo) G and (vo;co) G ^4, then (w(t);m(t)) G By 2 . Hence we derive from the last 
estimate that 



w(i) ■ Vm(t),*(t) 



<7- 1 (W2C(fi)) 2 |m(t)| 2 + 7 ||*(i)||?. 



(68) 



QUALITATIVE PROPERTIES OF THE MHD EQUATIOSN FOR BIPOLAR FLUIDS 



18 



Finally, 



8(t)- Vm(t),*(t) 



< C(n)||e(t)|[ 2 |*(t)| (HmCOIl! + |m(t)|) . 



Since (w(i);m(i)) £ -By 2 , we infer from the last estimate that 



6(t)-Vm(i),$(t) 



<^ 2 <7(0)||Q(i)|| 2 |*(i)|, 



(69) 



Using dSSJ-dSSD in jSSJ we deduce 
ld|^(t)| 2 



2 dt 



+ s\\*(t)\\ 2 < 2x||e(t)||l + 37||*(t)||? + (mC(ft)) 2 (iox- 1 + 7- 1 ) |*(t)| 2 

+ 2 7 - 1 ( m C(fi)) 2 |m(t)| 2 . (70) 

Letting Y(t) = |Q(i)| 2 + l*^)) 2 and summing up (f64"|) and (f70|) side by side implies that 
1 dY{t) 



2 eft 



+ ^\mm + s\\^{t)\\i < (mcwy \\w( t )u + 7x||G(t)n 2 + Hmm 

+ ( W2 C(ft)) 2 (Ux- 1 + |^)( 2 + x- 1 (C(ft)) 2 |Q(i)| 2 

( W2 C(fi)) 2 |w(i)| 2 + (fi P2 C(n)) 2 (x- 1 + 27- 1 ) |m(i)| 2 . 

(71) 

s 

Choosing x = — and 7 = — . We deduce from (|7ip that there exist positive constants 
14 6 

Ci,i = 1, ... ,5, depending only on U, [i, S, /Ui, p 2 such that 

Y^T + fH WH 2 + f H*(*)ll? ^ ^l WI 2 + &l*(*)l 2 + ^|w(i)| 2 + C 4 ||w(t)|| 2 

+(5 5 |m(t)| 2 . 

(72) 

Let us set xi = 2max(C*i, (7 2 ), X2 = 2max(C , 3, C4), and X3 = 2max(C < 5, 1). We infer from 
(J72D that 



eft 



where 



+ m ||G|| 2 + 511* ( t )\\( < X Y(t) + X2|^(t)r + Xs\\Z 



|Z(t)| 2 = |w(t)| 2 + |m(t)| 2 , 



\Z(t)\\l = \\w(t)\\ 2 + \\ m ^^ 



\Z \ 2 = |w | 2 + |m 



Dropping out the second term in the left hand side of (j72|) and invoking Gronwall's lemma 
yield 



Y(t) < X Ae Xlt / {\Z(s)\ 2 + \\Z(s)\\ 2 Y ) ds,Vt > 0, 
J 

where X4 = max (X2,X3)- From (|46p we have 

|Z(t)| 2 < iZope^ 1 *, 
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and 

v I \\Z(s)\§ds< m [ \Z(s)\ 2 ds, 
Jo Jo 

< \Z \ 2 (e^ - 1) , 



< l^nlV 1 *. 



Therefore 



Y(t)< X ^ 1+Xl)t \Z \ 2 (- + -). 

\m v o) 



And this estimate implies that for any 

|w(t)-U(t)| 2 + |m(t)-M(t)p 



0, 



|uo - v | + I b - c | 

for any t > as |uo — vo| — > and \bbo — Co| — > 0. This means that the semigroup 
§(i) : M — > H is uniformly differentiable wrt to the initial data. □ 

5. Bounds for the fractal df(A) and Hausdorff da (A) dimensions of A 

The fractal dimension of A is written as dAA). It is basically based on the number of 
closed balls of a fixed radius 5 needed to cover A (see, for instance, [331 Section 13.1.1]). We 
denote by N(A, 5) the minimum number of balls in such a cover. We recall the following 
definition for sake of precision. 

Definition 5.1 (See, for e.g., [33] ). If the closure of A is compact, then the fractal 
dimension of A is given by 

log N (A, 5) 
d / (^)=hmsup^ lQg(m) , 

where we allow the limit in the above equation to take the value +oo. 

Now we recall the definition of the Hausdorff dimension of A. It is based on the 
approximation of the d— dimensional volume of A by a covering of finite number of balls 
B(xi,r{) with radiie r% < 5. Let 

fj,(A,d,8) = mil "^rf : n < 5 and A C UiB(xi,ri) L 

i 

Definition 5.2. The Hausdorff dimension of A is given by 

d H (A) = inf {d : lim a(A, d, 5) = 0}. 
d>0 (5-s-O 

We have the following relation between fractal and Hausdorff dimensions. 
Theorem 5.3. d H (A) < d f (A). 

To calculate the bounds for the fractal dimension df(A) of A we will mainly follow the 



scheme in [331 Section 13.2]. For this purpose we let j$j = y^.J : « = 1,2, . . . )■ be an 
orthonormal basis of EI and P m be an orthogonal projection from H onto 

Span|<l>i : i = 1, 2, . . . ,m 
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For a solution (u; b) of (|2j) we denote by £(u; b) the mapping denned by 

dt \M(t)J L[U ^> \M(t).J 

Following the argument in [HJ chapter V] we denote by {A, ;i € N} the set of Lyapunov 
exponents of H and we set 

1 f* 

q m = hm sup - / -tr(£(u;b) o P m )(s)ds. 

t ^ 00 (uo;b )G^ 1 JO 

It is a standard fact (see, for instance, |41|, Chapter V]) that 

Ai H h A m < -q m . 

Therefore, it follows from [331 Theorem 13.16] (See also, [HJ Subsection V.3.4]) that 
Theorem 5.4. If —q m < 0, then m is the smallest positive integer such that df(A) < m. 
The value of £(u;b) at (4>i,^i) T is defined by £(u; b)(0j, ipi) T = (X;Y), where 

(X, 0) = - m - (r(u)5(&) - aAijki (u(t)) 

+ hi(t)-V&,^ + ^-Vu(t),^ • //|r, -Vu • b(t) • V ). 

and 

(y, = s(A^) -JUi- Vb(t),A + (u(t) ■ v^,A - Ui ■ Vu(t),^ 

for any (</>; € V. Therefore 

£(u;b)(&;^) T ,(&5^)) = (X,^) + (y,^), 

which is equivalent to 
£(u;b)(&;^) T , =-^(^^'^^) " ( r (u)£(<« " a A ijkl S 

-S\Vi>i\ 2 + (^•Vu(t),^] (^i- V0i,b(t) ) + ( Vb(t),V, 



• Vu(t),t/), 



Let us recall Equation (2.16)] 

I = I (r(u(t))|£-(^)| 2 - aA ijkl (u(t))\Sm 2 )dx >2evo [ , Jf ^* 1+a/2 dx 

Jn\ J Jn (e + \£ (u(t))\ 2 )^ < 



n (e + |£(u(t))| 2 ) J 
+ 2(l-aVo / |g( ^ )|2 , 2 dx 
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Let us set I\ (resp., I 2 ) the first (resp., second) in the above inequality. Both of I\ and I 2 
are positive, so we have 

\T ,\\ ^ ( 9£((f)i)\ 2 



r^b)^;^,^;^) ) >Ml( -5^'^^ ) +^|V^| 2 + 2(l-a)/ 2 



• Vu(t), 0i )-//(( Vi " V^i, b(t) ) + ( <f>i ■ Vb(t), ^ 

• Vu(t),^ 



Since < a < 1, we easily see that p = 2 — a 6 (1, 2]. We can check that 

(2(l-a))- 1 / 2 > ^( £ +|^(u(t))| 2 ) p/2 ^ " QV(&)|*dxV /P 
from which along with Korn's inequality we derive that 

p-2 

(2(l-a))- 1 / 2 > H ( e +|5(u(*))|Y /2 dx) " ^)[|^[|^. 
For < a < 1 we have that 

f (s+\£( y u(t))\ 2 ) p/2 dx<C(n,a,e)+ [\£(u(t))\ p dx. (73) 
As 1 < p < 2 we also see that 

\P/2 



/ (e+ |£(u(i))| 2 ) p/2 dx < C(n,a,e) +C(fi)( / |£ (u(i))[ 2 dx J 

\Jn / 

Since (u; b) is a solution of ([2]), u is an element of L 2 (0, oo;Vi). Therefore 

/ (e + \£ (u(t))\ 2 ) p/2 dx < C(n,a,e) + C(n,a), 
Jn 

for almost everywhere t > 0. Thanks to the fact that p — 2 > 0, we see that 

p-2 

(2(l-a))- 1 / 2 > (<7(n,a,e)) P X 
where C(0,a,e) = C(0,a, e) + C(f2, a). Hence 

- tr(£(u; b) o P m ) = - £ (c(u; b)*„ *,) > Ml f (^M, ^M) + 5 f)| V* 

i=l ^ ' i=l ^ ' i=l 

m m 



(74) 



i=i i=i 



P-2 

where iT(Sl) =2(1 -a) I C(p,a,e) ) * and 



F(u, b, <j> i} Vi) = • Vu(t), 0*) + /i (Vi ■ b(t)) + • Vb(t), V<) + W>< • VVi, u(t)) . 
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From Holder's inequality and Young's inequality we derive that 

Vu(t),&) < | u(t) | 1 0* 1 1 1 0i [ | iyi , 

<xr lc, ( n ) 2 H u (*)ll2+^l^l 2 ll^llfri. 
for any h\ > 0. Since \4>i\ 2 + (V'il 2 = 1 we obtain from the last estimate that 

(&-Vu(t),&) < ^C^IIu^lll + xill^H^. 
By a similar argument we have that 

M f • V&,b(i)) + • Vb(t), < 2/xC(fi)|^|||^|| H i||b(t)|| 1 



<7^ 1 (2 / uC(fi)) 2 ||b(t)||2+ 7 ||^||| 2 , 

for any 7 > 0. Also 

H(il>i ■ VVi,u(t)) < ^C(fi)||Vi||i||u(t)|| 2 |^|, 

<^ 2 - 1 ( / xC(0)) 2 ||u(t)||| + x 2 ||^|| 2 ? i J 

for any x 2 > 0. Thus 

F(u,b,^,^) < xiH^II^ + 7 ||^||^ 2 +x 2 ||^|Ihi +C(fi) 2 (xfi + x^V 2 ) ||u(t)||l 

+ 7 - 1 (2^)) 2 ||b(t)|| 2 . 
This last inequality and Korn's inequality enable us to state that 

mm m 

-tr(£(u;b) o P rn ) > ^Kip) \\<t>i\\ 2 H 2 +SJ2 \\A\\m + E \\</h\\w^p 

i=l i=l i=l 

m m 

-tE H^Hh 2 - ^ E 11^11^ - 7- 1 (^)) 2 iib(t)in 

i=l i=l 



-X2 £ H^lllf! " C(ft) 2 (xf 1 + X 2 -V 2 ) ||u(i)|||. 



1=1 



It was proved in [5] that there exists a positive constant such that for any a > 



TM* iifa=0 ' 

' |W " P " l^jlHIl! (*¥) II* otherwise, 

with 8' = 2 (f^g). Therefore 

-tr(£(u;b)oP m ) < f ^#(0) - J E ll&lllp - tE I^Hh 2 

^ ' 8=1 2=1 



i=l i=l 
2 ( „-l i „-1..2\ IL,C + M|2 „,-lAy( ,,^/nV^2||u/+Ml2 



- C(ny (x^ + x^V) ||u(*)||^ - 7- i (4/xC(0))^||b(t)||f. 
As (7 and 7 are arbitrary, we can choose them in such a way that 
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/ n x K(n)S' \^ , uiK(n) ml , 

that is a = — and 7 = . This choice implies that 

\K(n)(l-6')J 2 



-tr(£(u;b)oP m ) > 



5>d(n) \K(n)(l-5') 



-1 ra m 



8// 2 



i=l i=l 

c(nf (xr 1 + x 2 -y) || U (t)ii3 - _^||b(t)nf. 



Choosing x 2 = S and k\ = %-, where 



2 ' 

1 



7 S'd(n)\K(n)(i-S')J ' 1 ] 



we infer from the last estimate that 



1 m / 9 2 \ 

tr(£(u; b) o Pm ) > 2- £ H^ll^ - C(tt) 2 ( 5" 1 + ||u(i)|| 2 

8=1 \ 7 / 



8/i 



2 



|b(t)ll! 



Choosing as the eigenfunctions of the Stokes operator whose eigenvalues satify 

Xj > cXd n , (76) 

we derive from |33] that 

m 

Y,mm>^im 1+2/n . (77) 



i=l 

Thus 



-tr(£(u;b)oP m ) > 2^i m i + 2/n_ c(fi)2 f s -i + *£\ Mm 2_ 



' V "b(t)||l. (78) 



To give the bounds for the fractal dimension of A we also need to estimate the time average 
of \\u(t)\\% and ||b(i)||f. From (Q3D we derive that 

|u(i)| + |b(i)|+2/ii#(ft) f t \\u(s)\\ 2 2 ds + 2S / < ||b(s)||?^<47" 1 |/| 2 + 7 /"'|u(a)| 2 ek 

Jo Jo Jo 

+ |u | 2 + jb | 2 . 

By using Poincare's inequality and choosing 7 = ^ 1 — — -we can deduce from the above 

2Ai 

inequality that 

If 1 A 
lim - / ||u(s)|||ds < — — -, 

lim- /W)lli^<^, 
t-toa t Jq b 



where 



Now we see that 
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Owing to [33j Theorem 13.16] we see that df(A) < m if m is the smallest positive integer 
such that —q m < 0, that is 



m — 1 < 



2A 



c 2 (n) 



1 2/x 2 

+ ' 



+ 



8^ 



n + 2 



< m, 
2 -a 



(81) 



ic\i^K{^i) \ v ' \s i j s 

where 7' and A are given respectively by (1751) and (1791) . and 5' = 2 

\4 + oCj 

The conclusion of this section is stated in the following claim that we have just proved 
above. 

Theorem 5.5. The global attractor A for the MHD equations for the nonlinear bipolar 
fluids (|2|) is finite- dimensional with da (A) < df(A) < m where m > is positive integer 
given by 



m — 1 < 



2A 



1 2fj, 2 \ 8fi 



n + 2 



ic\imK{VL) \ x ' \s i ) s 

where 7' and A are given respectively by ([75]) and ([79]) . ana" 5' = 2 



< m, 
2 - a 



(82) 



4 + a 
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